Desynchronization waves and localized instabilities in oscillator arrays 
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We consider a ring of identical or near identical coupled periodic oscillators in which the con- 
nections have randomly heterogeneous strength. We use the master stability function method to 
determine the possible patterns at the desynchronization transition that occurs as the coupling 
strengths are increased. We demonstrate Anderson localization of the modes of instability, and 
show that such localized instability generates waves of desynchronization that spread to the whole 
array. Similar results should apply to other networks with regular topology and heterogeneous 
connection strengths. 
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In this Letter we discuss the synchronization of a large 
number of near-identical oscillators that are locally cou- 
pled with connections of random strength. Synchroniza- 
tion in networks of coupled oscillators has recently re- 
ceived considerable interest 1] , and has relevance in fields 
like biology , chemistry (sf , lasers P, Q , and commu- 
nications 13 . Usually, the networks studied have been 
assumed to have connections of equal strength. In prac- 
tice, the connections between different oscillators may 
have different strengths, and in some cases this strength 
could have a large spread (e.g., in biological systems). A 
model and analysis method has been proposed by Pec- 
ora and Carroll XJ to systematically determine the sta- 
bility of the synchronized state in a network of identi- 
cal coupled oscillators. This method, the master stabil- 
ity function, has been used to study the synchronization 
properties of different networks @, ^ . Deng et al. [lol | 
have obtained, using the master stability function tech- 
nique, conditions for the distribution of the connection 
strengths that yield average stability of the synchronized 
state. Galias and Ogorzalek 11] have studied the effect 
of adding small perturbations to the coupling strengths 
in relatively small arrays of coupled chaotic oscillators. 
Denker et. al. J^] have studied the effect of small cou- 
pling strength heterogeneity in networks of pulse-coupled 
oscillators. Our approach in this Letter will be different: 
we consider the coupling strengths to have a relatively 
large spread, and will discuss phenomena that can be ex- 
pected when a large number of periodic oscillators are 
coupled in such a network. In particular, we will see 
that as the coupling strength is increased, the oscillators 
desynchronize in a localized region. The localization re- 
sults because the connection matrix has random compo- 
nents and the eigenvectors of this matrix are Anderson 



localized [isL IT^ . The effect of the localized instabil- 
ity spreads as a wave throughout the array, eventually 
resulting in an ordered state. Remarkably, in the case 
where the oscillators are not identical the final state of 
the locally unstable system is more ordered than in the 
case where the system is stable. 

We consider a model system of N identical dynamical 
units, each one of which, when isolated, satisfies Xi = 
F{Xi), where i = 1,2, ... N, and Xi is the d-dimensional 
state vector for unit i. (The case of nearly identical units 
is considered at the end of this Letter. See also 
The oscillators, when coupled, are taken to satisfy (e.g.. 



N 



X,^F(X,)-gY,G^JH{X, 



(1) 



where the coupling function H is independent of i and 
j, and the matrix G is a symmetric Laplacian matrix 
Gij = 0) describing the network connections. The 
constant g determines the global strength of the coupling. 

There is an exactly synchronized solution of Eqs. 
Xi ~ X2 ^ ■ ■ ■ ~ Xn = s{t), whose time evolution is 
the same as the uncoupled dynamics of a single unit, s = 
F{s). In this Letter we will be concerned with the case 
where the synchronized state is periodic, s{t + T) = s{t). 
The stability of the synchronized state can be determined 
from the variational equations obtained by considering an 
infinitesimal perturbation from the synchronous state, 
X,(t) = s(t) + e,(t). 
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Let e = [ei, 62, . . . , Cat], and define the d x N ma- 
trix T] = [771, 772, . . . , 77Ar] by e = i']L^ , where L is the 
orthogonal matrix whose columns are the correspond- 
ing real orthonormal eigenvectors of G; GL = LA, 
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A = (iiag(Ai, A2, . . . , Aat) where is the eigenvalue of 
G for eigenvector k. Then Eqs. are equivalent to 



m - {DF{s)- gXkDH{s))T^k. 



(3) 



The quantity rjk is the weight of the fc*'* eigenvector of G 
in the perturbation e. The linear stability of each 'spatial' 
mode k is determined by the stability of the zero solution 
of By introducing a scalar variable a = gXk, the 
set of equations given by lO can be encapsulated in the 
single equation, 



iDF{s) - aDH{s))ri. 



(4) 



The master stability function ^'(a) is the largest Lya- 
punov exponent for this equation. This function de- 
pends only on the coupling function H and the chaotic 
dynamics of an individual uncoupled element, but not 
on the network connectivity. The network connectiv- 
ity determines the eigenvalues Afe (independent of details 
of the dynamics of the chaotic units). The stability of 
the synchronized state of the network is determined by 
= sup;, '^{gXk), where > indicates instability. 
As an illustrative example, we consider periodic 
Rossler oscillators, obeying the equations 



y = x + 0.2?/, 
z = 0.2 + z{x -2.5). 



(5) 



In terms of our previous notation, d = 3, and X = 
[x, y, z]^. The master stability function for this system is 
shown in Fig. ^ As seen in this figure, 5* (a) approaches 
zero from negative values as a — > 0+. This is a general 
feature for systems where the individual, uncoupled units 
are stable limit cycle oscillators. We also see that 'i'{a) 
crosses from negative (stable) values to positive (unsta- 
ble) values at a critical a value {a w 4.15). The existence 
of such a transition is a robust feature that depends on 
the type of coupling and oscillator. We now consider a 




FIG. 1: Master stability function "if (a) versus a for Eqs. 
At Q ~ 4.15, the master stability function becomes positive. 



network of N of these oscillators nearest-neighbor cou- 
pled in a ring, such that the strength of each individ- 
ual link is random. The coupling strengths are obtained 
from an independent and identically distributed random 
sequence {ai}fLi. The matrix G is then 



G 



/ bi -ai ••• ~aN \ 
-ai 62 -02 ••• 

—02 ^3 —03 • • • 



\ —aN — ajv-i 



(6) 



where bi = (aj_i-t-aj) for i = 1, . . . , N (we take og = a^y). 

The eigenvectors of the matrix G determine the pos- 
sible desynchronization patterns. It is known that 
the eigenvectors of certain types of random matrices 
are exponentially localized (e.g., Anderson localization 
USQl)- In our case, the eigenvector {ui}fLi with eigen- 
value A satisfies 



ti+i = a.j+i(A -I- flj -I- ai+i - Qit.-^ -^), 



(7) 



Viewing Eq. Q as a random dynamical 



where i,- = - 

system for ti, we find numerically that in our case, 



1 " 

lim -y^log(|tj 



(8) 



exists and is independent of the initial condition and 
noise realization. Eigenvectors of © tend to have a lo- 
calized amplitude peak at some location zq and decay like 
\ui\ oc e'''l*~'°l away from the peak 7"-'^ is thus the 
localization length. 

We choose the a^'s to be uniformly distributed in 
(0.1, 1) (note that any multiple of this would lead to 
the same eigenvectors). (Since ai > 0.1 we avoid the 
possibility ^ 1 that would effectively disconnect the 
network.) The effects we will describe for this network 
should be regarded as an example of what could be ex- 
pected in more general networks with random coupling. 
In Fig. |2Ia) we show the eigenvector with largest eigen- 
value for a realization of the matrix G using N = 500. 
Figure |2tb) shows the localization length 7"^ as a func- 
tion of A calculated using Eq. (jS)). The eigenvectors are 
seen to be sharply localized for the largest eigenvalues, 
and become less localized as the eigenvalues decrease. 

As the coupling strength g is increased, the eigenvec- 
tors with largest eigenvalue become unstable. These 
eigenvectors have the smallest localization length (see 
Fig. 121 (b)). We will now describe what occurs in this 
situation. We fixed the same realization of the matrix G 
used in producing Fig. |2Ja) . The four largest eigenvalues 
are 3.61, 3.41, 3.38, and 3.30. For g = 1.24 the eigen- 
vector with largest eigenvalue is unstable, and the next 
two eigenvectors are barely unstable {a = 4.47, 4.23 and 
4.19 in Fig. (^)). We start with initial conditions near 
the synchronized state and then let the system evolve 
according to Eqs. (^). In Fig. |2|we show snapshots of Xi 
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FIG. 2: (a) Eigenvector Ui for the largest eigenvalue A = 
3.61 for a particular realization of the matrix G in with 
A'' — 500. It is sharply localized. Note that the range is only 
200 < i < 300. (b) Localization length 7"^ calculated using 
Eq. ©. 




and n{i,t) is an integer chosen so that </) is a continuous 
function of t and that + 1, t) is close to t) for all 
i. Figure^ shows two snapshots of the x coordinate and 
the phase as defined above as a function of i (the i origin 
was displaced so that what happens opposite the location 
of the unstable mode can be observed clearly, and for 
each time a constant was added to (/> so that max^ (/) — 
0). As can be observed in the Figs. ^ (a) and (c), a 
region with a constant phase gradient expands on both 
sides of the unstable mode. In the final state (Figs.^l^b) 
and (d)) the phase has a minimum at the location of 
the unstable mode and increases linearly on both sides 
reaching a maximum at the opposite end of the ring. 
This phase profile increases uniformly with time. The 
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FIG. 4: Plots (a) and (b) shows the x coordinate of oscillator 
i versus the site index i for times 3750 and 9660. Plots (c) 
and (d) show the phase of oscillator i at the same times as 
for (a) and (b) respectively. Compare with Eq. (llOt 



FIG. 3: Plots of the x coordinate of oscillator i versus the site 
index i, at times (a) 0, (b) 1400, (c) 2800, (d) 4200, (e) 5600, 
and (f) 10000. All the plots have the same scale as (e). 



as a function of the site index i for six successively in- 
creasing times. Starting from a nearly synchronized state 
(Fig. O^a)), the oscillators desynchronize at the location 
(see Fig. |2;a)) of the localized mode (Fig. Elb)). The 
desynchronization spreads as a wave to farther regions 
of the array (Figs. |2tc)-(e)). At the end, the domain of 
the wave covers the entire array (Fig. |2If)). This pro- 
cess is dominated by the most unstable mode. The other 
two less unstable modes can be seen as tiny defects at 
i « 327, 402 in the otherwise smooth wave. (The effect 
of these less unstable modes is most evident in Fig.l^Jc). 
They also have a discernible, although small, effect in the 
final state [arrows in Fig. |2If)].) 

The final state and the process leading to it can be 
understood in terms of the phase of the oscillators. Define 
the phase </>(«, t) = 27r{n(i,t) + {t — t_(i, t))(t+(i, t) — 
t_(i,t))^^}, where t^{i,t) = max{s : Xi{s) = 0,ii > 
0, s < t}, t+{i,t) = min{s : Xi{s) = 0, ii > 0,s > t}. 



cause of this phenomenon is that, as the oscillators in 
the region of the unstable mode desynchronize, they go 
to limit cycles that have a slightly lower frequency than 
that of the original orbit. Oscillating at a slower pace 
than the others, they drag the adjacent oscillators, and 
these drag theirs in turn, continuing until an equilibrium 
is reached. An equation describing approximately the 
evolution of the phase of the oscillator at location ^ and 
time i, 0(^, t), in a chain of diffusively coupled oscillators 
is given in the continuous limit by \w{ 



dt 



MO, 



(9) 



where w(^) is the frequency of the oscillator at location 
^, and a and b are constants. If this frequency is suffi- 
ciently smaller (larger) in a localized region and b is neg- 
ative (positive) , the equation predicts the development of 
waves that emanate from that region. The phase profile 
resulting from such forcing in a small region centered at 
the orig in (|^| < I) can be approximated for large ^ and 

t as ITgi 



0(e,i) =w;oi-max(0,/c(t;i-|e|)), 



(10) 
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where wq = w{^) for |^| > / and k and v depend on a and 
b and w{^). For appropriate k and w, equation H1U|I agrees 
well with Figs.^(c) and (d). In the example presented 
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FIG. 5: Each plot shows the x coordinate of oscillator i as 
a function of the site index i. The time is 0, 1400, 2800, 
4200, 5600, and 9970 for plots (a) to (f) and similarly for 
plots (g) to (1). A parameter mismatch was introduced in the 
oscillators, (a)-(f): All the modes are stable, (g)-(l): The 
pattern is organized by an unstable mode as in Fig. |HJf ) . All 
the plots have the same scale as (e). 

above, the pattern created by the unstable mode can be 
regarded as a more disordered synchronization than that 
of the original identical synchronization. However, in re- 
alistic situations, an unstable mode can actually make 
synchronization more orderly. In real systems, small dif- 
ferences in the parameters or small noise are expected. 
Under these circumstances, the different oscillators will 
be subject to small perturbations. The modes with eigen- 
values close to zero have a master stability function close 



to zero (see Fig. and also are nearly unlocalized (see 
Fig.EIb)). Thus, the phase of each oscillator will be sub- 
ject to perturbations whose projection onto the nearly 
unlocalized modes are only very weakly damped. The 
identical synchronization of the array is thus spoiled by 
mismatch or noise. As an illustration, we randomly per- 
turb the parameters of the different oscillators, so that 
they lie within ±3% of the original parameters. We then 
solved Eqs. Q with g — 1.1 and g — 1.24. For g = 1.1, 
all the modes are stable; in the case g = 1.24, three 
modes are stable as discussed above. In Figs. El (a)- (f) 
we show snapshots of the case g — 1.1, and in Figs. 
(g)-(l) we show the corresponding snapshots for the case 
g — 1.24. When all of the modes are stable, the system 
exhibits a state in which there is erratic slow variation 
of the Xi with i. When there is an unstable mode, how- 
ever, a more organized state is reached. If one picks two 
different oscillators j and fc, they will satisfy asymptoti- 
cally Xj (t — r) = Xk (t) , where r is a simple function of j 
and k (see Fig. 111(d) ). Thus the oscillators are pairwise 
lag synchronized |17l | . In realistic large arrays of periodic 
oscillators, it might be convenient to have one unstable 
mode. This mode will, despite its localized nature, in- 
duce global organization of the system (Fig. 01. 

In conclusion, we find that large arrays of periodic os- 
cillators locally coupled by connections of randomly het- 
erogeneous strength can experience a desynchronization 
transition characterized by the appearance of unstable 
Anderson localized modes. Furthermore, we find that, 
past the transition, the localized mode plays the key role 
in organizing the final global pattern of the system oscil- 
lations. 
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